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ABSTRACT: 

Magnetic fields in galaxy clusters have been measured using a variety of techniques, including: 
studies of synchrotron relic and halo radio sources within clusters, studies of inverse Compton 
X-ray emission from clusters, surveys of Faraday rotation measures of polarized radio sources 
both within and behind clusters, and studies of Cluster Cold Fronts in X-ray images. These 
measurements imply that most cluster atmospheres are substantially magnetized, with typical 
field strengths of order 1 /iGauss with high areal filling factors out to Mpc radii. There is likely, 
however, to be considerable variation in field strengths and topologies both within and between 
clusters, especially when comparing dynamically relaxed clusters to those that have recently 
undergone a merger. In some locations, such as the cores of cooling flow clusters, the magnetic 
fields reach levels of 10-40 fiG and may be dynamically important. In all clusters the magnetic 
fields have a significant effect on energy transport in the intracluster medium. We also review 
current theories on the origin of cluster magnetic fields. 

1 Introduction 

Magnetic fields play an important role in virtually all astrophysical phenomena. 
Close to home, the Earth has a bipolar magnetic field with a strength of 0.3 
G at the equator and 0.6 G at the poles. This field is thought to originate in a 
dynamo due to fluid motions within the liquid core (Soward 1983). With its faster 
angular rotation, Jupiter leads the planets with an equatorial field strength of ~ 
4 G (Warwick 1963; Smith et al. 1974). A similar mechanism produces the solar 
magnetic fields which give rise to spectacular sunspots, arches, and flares (Parker 
1979). Within the interstellar medium, magnetic fields are thought to regulate 
star formation via the ambipolar diffusion mechanism (Spitzer 1978). Our own 
Galaxy has a typical interstellar magnetic field strength of ~2 jiG in a regular 
ordered component on kiloparsec scales, and a similar value in a smaller-scale, 
random component (Beck et al. 1996; Kulsrud 1999). Other spiral galaxies have 
been estimated to have magnetic field strengths of 5 to 10 [iG, with fields strengths 
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up to 50/iG found in starburst galaxy nuclei (Beck et al. 1996). Magnetic fields are 
fundamental to the observed properties of jets and lobes in radio galaxies (Bridle 
& Perley 1984), and may be primary elements in the generation of relativistic 
outflows from accreting, massive black holes (Begelman, Blandford, & Rees 1984). 
Assuming equipartition conditions apply, magnetic field strengths range from a 
few /^G in kpc-scale extended radio lobes, to mG in pc-scale jets. 

The newest area of study of cosmic magnetic fields is on larger scales still, that 
of clusters of galaxies. Galaxy clusters are the largest virialized structures in the 
universe. The first spatially resolving X-ray observations of clusters (Forman et 
al. 1972) revealed atmospheres of hot gas (10 7 to 10 8 K) which extend to Mpc 
radii and which dominate the baryonic mass of the systems (10 13 to 10 14 Mg). 
Soon thereafter came the first attempts to measure magnetic field strengths in the 
intracluster medium (ICM) (Jaffe 1977). Only in the last decade has it become 
clear that magnetic fields are ubiquitous in cluster atmospheres, certainly playing 
a critical role in determining the energy balance in cluster gas through their 
effect on heat conduction, and in some cases perhaps even becoming important 
dynamically. 

Cluster magnetic fields have been treated as secondary topics in reviews of 
cluster atmospheres (Sarazin 1988; Fabian 1994), and in general reviews of cosmic 
magnetic fields (Kronberg 1996; Ruzmaikin, Shukurov, & Sokolov 1987). To date 
there has been no dedicated review on cluster magnetic fields. 

The focus of this review is primarily observational. We summarize and critique 
various methods used for measuring cluster magnetic fields. In the course of the 
review we consider important effects of magnetic fields in clusters, such as their 
effect on heat conduction and gas dynamics, and other issues such as the lifetimes 
of relativistic particles in the ICM. We then attempt to synthesize the various 
measurements and develop a general picture for cluster magnetic fields, with the 
caveat that there may be significant differences between clusters, and even within 
a given cluster atmosphere. We conclude with a section on the possible origin of 
cluster magnetic fields. 

We assume Hg = 75 km s _1 Mpc -1 and qn=0.5, unless stated otherwise. Spec- 
tral index, a, is defined as S v oc v a . 

2 Synchrotron radiation 

2.1 Radio halos 

Over 40 years ago Large (1959) discovered a radio source in the Coma cluster 
that was extended even when observed with a 45' beam. This source (Coma 
C) was studied by Willson (1970) who found that it had a steep spectral index 
and could not be made up of discrete sources, but instead was a smooth "radio 
halo" with no structure on scales less than 30'. Willson further surmised that 
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the emission mechanism was likely to be synchrotron, and if in equipartition 
required a magnetic field strength of 2 /xG. In Fig. 1 we show the best image 
yet obtained of the radio halo in the Coma cluster. Other radio halos were 
subsequently discovered, although the number known remained under a dozen 
until the mid-90s (Hanisch 1982). 

Using the Northern VLA Sky Survey (NVSS; Condon et al. (1998)) and X-ray 
selected samples as starting points Giovannini & Feretti (2000); Giovannini, Tordi 
& Feretti (1999) have performed moderately deep VLA observations (integrations 
of a few hours) which have more than doubled the number of known radio halo 
sources. Several new radio halos have also been identified from the Westerbork 
Northern Sky Survey (Kempner & Sarazin 2001). These radio halos typically 
have sizes ~1 Mpc, steep spectral indices (a < —1), low fractional polarizations 
(< 5%), low surface brightnesses (~ 10 -6 Jy arcsec -2 at 1.4 GHz), and centroids 
close to the cluster center defined by the X-ray emission. 

A steep correlation between cluster X-ray and radio halo luminosity has been 
found, as well as a correlation between radio and X-ray surface brightnesses in 
clusters (Liang et al. 2000; Feretti et al. 2001; Govoni et al. 2001a). A complete 
(flux limited) sample of X-ray clusters shows only 5% to 9% of the sources are 
detected at the surface brightness limits of the NVSS of 2.3 mjy beam -1 , where 
the beam has FWHM = 45" (Giovannini & Feretti 2000; Feretti et al. 2001). 
But this sample contains mostly clusters with X-ray luminosities < 10 45 erg 
s . If one selects for clusters with X-ray luminosities > 10 45 erg s _1 , the radio 
detection rate increases to 35% (Feretti et al. 2001; Owen, Morrison & Vogues 
1999). Likewise, there may be a correlation between the existence of a cluster 
radio halo and the existence of substructure in X-ray images of the hot cluster 
atmosphere, indicative of merging clusters, and a corresponding anti-correlation 
between cluster radio halos and clusters with relaxed morphologies, e.g., cooling 
flows (Govoni et al. 2001a), although these correlations are just beginning to be 
quantified (Buote 2001). 

Magnetic fields in cluster radio halos can be derived assuming a minimum 
energy configuration for the summed energy in relativistic particles and magnetic 
fields (Burbidge 1959), corresponding roughly to energy equipartition between 
fields and particles. The equations for deriving minimum energy fields from radio 
observations are given in Miley (1980). Estimates for minimum energy magnetic 
field strengths in cluster halos range from 0.1 to 1 /xG (Feretti 1999). One of the 
best studied halos is that in Coma, for which Giovannini et al. (1993) report a 
minimum energy magnetic field of 0.4 /j,G. These calculations typically assume 
k = 1, r] = 1, fi ow = 10 MHz, and !4,; g h = 10 GHz, where k is the ratio of 
energy densities in relativistic protons to that in electrons, r\ is the volume filling- 
factor, i\ ow is the low frequency cut-off for the integral, and z^hjgh is the high 
frequency cut-off. All of these parameters are poorly constrained, although the 
magnetic field strength only behaves as these parameters raised to the S power. 
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Figure 1 : WSRT radio image of the Coma cluster region at 90 cm, with angular resolution 
of 55" x 125" (HPBW, RA x DEC) from Feretti et al (1998). Labels refer to the halo source 
Coma C and the relic source 1253+275. The grey scale range displays total intensity emission 
from 2 to 30 mjy/beam while contour levels are drawn at 3, 5, 10, 30, and 50 mjy/beam. The 
bridge of radio emission connecting Coma C to 1253+275 is resolved and visible only as a region 
with an apparent higher positive noise. The Coma cluster is at a redshift of 0.023, such that 
1' = 27 kpc. 
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For example, using a value of k ~ 50, as observed for Galactic cosmic rays (Meyer 
1969), increases the fields by a factor of three. 

Brunetti et al. (2001a) present a method for estimating magnetic fields in the 
Coma cluster radio halo independent of minimum energy assumptions. They 
base their analysis on considerations of the observed radio and X-ray spectra, 
the electron inverse Compton and synchrotron radiative lifetimes, and reason- 
able mechanisms for particle reacceleration. They conclude that the fields vary 
smoothly from 2±1 in the cluster center, to 0.3±0.1 fiG at 1 Mpc radius. 

2.2 Radio relics 

A possibly related phenomena to radio halos is a class of sources found in the 
outskirts of clusters known as radio relics. Like the radio halos, these are very 
extended sources without an identifiable host galaxy (Fig. 1). Unlike radio halos, 
radio relics are often elongated or irregular in shape, are located at the cluster 
periphery (by definition), and are strongly polarized, up to 50% in the case of 
the relic 0917+75 (Harris et al. 1993). As the name implies, one of the first 
explanations put forth to explain these objects was that these are the remnants 
of a radio jet associated with an active galactic nucleus (AGN) that has since 
turned off and moved on. A problem with this model is that, once the energy 
source is removed, the radio source is expected to fade on a timescale << 10 8 
years due to adiabatic expansion, inverse Compton, and synchrotron losses (see 
§4.1). This short timescale precludes significant motion of the host galaxy from 
the vicinity of the radio source. 

A more compelling explanation is that the relics are the result of first or- 
der Fermi acceleration (Fermi I) of relativistic particles in shocks produced dur- 
ing cluster mergers (Ensslin et al. 1998), or are fossil radio sources revived by 
compression associated with cluster mergers (Ensslin & Gopal-Krishna 2001). 
Equipartition field strengths for relics range from 0.4 - 2.7 hfJJ uG (Ensslin et 
al. 1998). If the relics are produced by shocks or compression during a cluster 
merger, then Ensslin et al. (1998) calculate a pre-shock cluster magnetic field 
strength in the range 0.2-0.5 fiG. 

3 Faraday rotation 

3. 1 Cluster center sources 

The presence of a magnetic field in an ionized plasma sets a preferential direction 
for the gyration of electrons, leading to a difference in the index of refraction for 
left versus right circularly polarized radiation. Linearly polarized light propa- 
gating through a magnetized plasma experiences a phase shift of the left versus 
right circularly polarized components of the wavefront, leading to a rotation of 
the plane of polarization, Ax = RM A 2 , where Ax is the change in the position 
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angle of polarization, A is the wavelength of the radiation, and RM is the Fara- 
day rotation measure. The RM is related to the electron density, n e , and the 
magnetic field, B, as: 

L 

RM = 812 J n e B • d\ radians m -2 , (1) 

o 

where B is measured in /iGauss, n e in cm -3 and d\ in kpc, and the bold face 
symbols represent the vector product between the magnetic field and the direc- 
tion of propagation. This phenomenon can also be understood qualitatively by 
considering the forces on the electrons. 

Synchrotron radiation from cosmic radio sources is well known to be linearly po- 
larized, with fractional polarizations up to 70% in some cases (Pacholczyk 1970). 
Rotation measures can be derived from multifrequency polarimetric observations 
of these sources by measuring the position angle of the polarized radiation as a 
function of frequency. The RM values can then be combined with measurements 
of n e to estimate the magnetic fields. Due to the vector product in Eq. 1, only 
the magnetic field component along the line-of-sight is measured, so the results 
depend on the assumed magnetic field topology. 

Most extragalactic radio sources exhibit Faraday rotation measures (RMs) of 
the order of 10's of rad m -2 due to propagation of the emission through the 
interstellar medium of our galaxy (Simard-Normandin, Kronberg & Button 1981). 
Sources at Galactic latitudes < 5° can exhibit ~300 rad m -2 . For the past 30 
years, however, a small number of extragalactic sources were known to have far 
higher RMs than could be readily explained as Galactic in origin. Large intrinsic 
RMs were suspected, but the mechanism(s) producing them were unclear. 

Mitton (1971) discovered that the powerful radio galaxy Cygnus A had large, 
and very different RMs (35 vs —1350 rad m -2 ), in its two lobes. While its low 
galactic latitude (5.8°) could possibly be invoked to explain the high RMs, the 
large difference in RMs over just 2' was difficult to reproduce in the context of 
Galactic models (Alexander, Brown, & Scott 1984). This "RM anomaly" was 
clarified when Dreher, Carilli & Perley (1987) performed the first high resolution 
RM studies with the VLA and found complex structure in the RM distribution 
on arcsec scales (Fig. 2), with gradients as large as 600 rad m -2 arcsec -1 . These 
large gradients conclusively ruled out a Galactic origin for the large RMs. 

Perhaps just as important as the observed RM structure across the lobes of 
Cygnus A was the discovery that the observed position angles behave quadrati- 
cally with wavelength to within very small errors over a wide range in wavelengths 
(Dreher, Carilli & Perley 1987). Examples of this phenomenon are shown in 
(Fig. 5). Moreover, the change in position angle from short to long wavelengths 
is much larger than n radians in many cases, while the fractional polarization 
remains constant. This result is critical for interpreting the large RMs for cluster 
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Figure 2: The RM distribution in Cygnus A based on multi-frequency, multi-configuration 
VLA observations. The resolution is 0.35" (Dreher et al. 1987). The colorbar indicates the 
range in RMs from —3400 to +4300 rad m -2 . Note the undulations in RM on scales of 10-30 
kpc. Contours are overlaid from a 5 GHz total intensity image. The RM was solved for by 
fitting for the change in polarization angle with frequency on a pixel-by-pixel basis (see Fig. 5). 



center radio sources, providing rigorous proof that the large RMs cannot be due 
to thermal gas mixed with the radio emitting plasma (Dreher, Carilli & Perley 
1987). Such mixing would lead to rapid depolarization with increasing wave- 
length, and departures from a quadratic behavior of \ with wavelength (Burn 
1966). 

The Cygnus A observations were the first to show that the large RMs must 
arise in an external screen of magnetized, ionized plasma, but cannot be Galactic 
in origin. Dreher, Carilli & Perley (1987) considered a number of options for the 
Faraday screen toward Cygnus A, and concluded that the most likely site was 
the X-ray emitting cluster atmosphere enveloping the radio source (Fabbiano et 
al. 1979). They found that magnetic fields in the cluster gas of 2-10 \xG could 
produce the observed RMs. 

Since the ground-breaking observations of Cygnus A, RM studies of cluster 
center radio sources have become a standard tool for measuring cluster fields. 
RM studies of radio galaxies in clusters can be divided into studies of cooling-flow 
and non-cooling-flow clusters. Cooling-flow clusters are those in which the X-ray 
emission is strongly peaked at the center, leading to high densities and cooling 
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times of the hot ICM in the inner ~100 kpc of much less than the Hubble time. 
To maintain hydrostatic equilibrium, an inward flow may be required (Fabian, 
Nulsen, & Canizares 1991). Typical mass cooling flow rates are 100 Mq yr _1 . 
The actual presence of material "cooling" and "flowing" is a topic that is hotly 
debated at present (Binney 2001). What is more agreed upon is that cooling-flow 
clusters are more dynamically relaxed than non-cooling flow clusters which often 
show evidence of cluster mergers (Markovic & Eilek 2001). 

Radio galaxies in cooling flow clusters attracted some of the first detailed RM 
studies by virtue of their anomalously high RMs (e.g., A1795: (Ge & Owen 1993); 
Hydra A: (Taylor & Perley 1993)). Out of a sample of 14 cooling-flow clusters 
with strong embedded radio sources Taylor, Barton & Ge (1994); Taylor et al. 
(2001b), found that 10 of 14 sources display RMs in excess of 800 rad m -2 , two 
(PKS 0745-191 and 3C84 in Abell 426) could not be measured due to a lack of 
polarized flux, and two (Abell 119 (Feretti et al. 1999) and 3C 129(Taylor et al. 
2001a)) have lower RMs, but with better X-ray observations turn out not to be 
in cooling-flow clusters. Hence, current data are consistent with all radio galaxies 
at the center of cooling flow clusters having extreme RMs, with the magnitude 
of the RMs roughly proportional to the cooling flow rate (see Fig. 3). 

The RM distributions for radio sources found at the centers of cooling flow 
clusters tend to be patchy with coherence lengths of 5-10 kpc (Fig. 4). Larger 
"patches" up to 30 kpc are seen for example in Cygnus A (Fig. 2). In both 
Cygnus A and Hydra A one can find "bands" of alternating high and low RM 
(see Figures 2 and 4). Such bands are also found in the non-cooling flow cluster 
sources (Eilek & Owen 2001), along with slightly larger coherence lengths of 15- 
30 kpc. In Hydra A there is a strong trend for all the RMs to the north of the 
nucleus to be positive and to the south negative. To explain this requires a field 
reversal and implies a large-scale (100 kpc) ordered component to the cluster 
magnetic fields in Hydra A. Taylor & Perley (1993) found the large scale field 
strength to be ~7 /iG, and more tangled fields to have a strength of ~40 /iG. 
A similar RM sign reversal across the nucleus is seen in A1795, although in this 
case the radio source is only 11 kpc in extent. 

Minimum cluster magnetic field strengths can be estimated by assuming a 
constant magnetic field along the line-of-sight through the cluster. Such estimates 
usually lead to magnetic field strengths of 5 to 10 /iG in cooling flow clusters, 
and a bit less (factor ~2) in the non-cooling flow clusters. 

From the patchiness of the RM distributions it is clear that the magnetic 
fields are not regularly ordered on cluster (Mpc) scales, but have coherence scales 
between 5 and 10 kpc. Beyond measuring a mean line-of-site field, the next level 
of RM modeling entails 'cells' of constant size and magnetic field strength, but 
random magnetic field direction, uniformly filling the cluster. The RM produced 
by such a screen will be built up in a random walk fashion and will thus have an 
average value of rad m , but a dispersion in the RM, orm, that is proportional 
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Figure 3: The maximum absolute RM plotted as a function of the estimated cooling flow rate, 
X, for a sample of X-ray luminous clusters with measured RMs from Taylor et al. (2001b). 
Both RM and X are expected to depend on density to a positive power, so in that sense the 
correlation is expected. 



to the square root of the number of cells along the line-of-sight through the 
cluster. The most commonly-fit form to the X-ray observations to obtain the 
radial electron density distribution, n e (r) through a cluster is the modified King 
model (Cavaliere & Fusco-Femiano 1976): 

n e (r)=n (l+r 2 /r 2 c )-W 2 , (2) 
where no is the central density, r c is the core radius, and (3 is a free parameter 



in the fit. Typical values for these parameters are 



200 kpc, /3 
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Figure 4: The RM distribution in Hydra A at a resolution of 0.3" (Taylor & Perley 1993) with 
total intensity contours overlaid. Multi-configuration VLA observations were taken at 4 widely 
spaced frequencies within the 8.4 GHz band, and a single frequency in the 15 GHz band. The 
colorbar indicates the range in RMs from —12000 to +5000 rad m~ 2 . 



For this density profile and cells of constant magnetic strength but random 
orientation, Felten (1996); Feretti et al. (1995) derived the following relation for 



Carilli & Taylor 



11 



-11016 




5 10 

CM**2 



5 10 
CM**2 



09 15 40.44 



40.43 



40.42 
RIGHT ASCENSION 



40.41 



40.40 



Figure 5: The observed position angles, x, of the linearly polarized radio emission as a function 
of the square of the observing wavelength, A 2 , for a number of positions in the southern lobe of 
Hydra A at a resolution of 0.3" (Taylor & Perley 1993). The points plotted are each separated 
by approximately one beamwidth and thus are independent of each other. This illustrates the 
consistency of the RMs within a coherence length of ~7 kpc. Notice also the excellent agreement 
to a A 2 -law for Ax = 600 degrees, nearly two complete turns. 



the RM dispersion: 

KB n r y V/ 2 / r(3/3-0.5) 

CT RM- (1+r2/r 2 )(6/3 _ 1)/4 y rm (6) 

where I is the cell size, r is the distance of the radio source from the cluster center, 
r is the Gamma function, and K is a factor which depends on the location 
of the radio source along the line-of-sight through the cluster: K = 624 if the 
source is beyond the cluster, and K = 441 if the source is halfway through the 
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cluster. Note that Eq. 3 assumes that the magnetic field strength, B, is related 
to the component along the line of sight, by B = y/3B\\. The cell size, 

I, can be estimated to first order from the observed coherence lengths of the 
RM distributions. Both cooling-flow and non-cooling flow clusters yield typical 
estimates of 5 to 10 kpc. Magnetic field strength estimates, however are two 
to three times higher in the cooling-flow clusters - 19 /^G in the 3C 295 cluster 
(Allen et al. 2001a) compared to 6 ^G in the 3C 129 cluster (Taylor et al. 2001a) 
using the methodology described above. 

Most radio sources found embedded in clusters are located at the center and 
identified with a cD galaxy. This relatively high pressure environment has been 
found in many cases to confine or distort the radio galaxy (Taylor, Barton & Ge 
1994), as well as giving rise to extreme RMs. For this same reason the extended 
radio sources in Hydra A and Cygnus A are unique in that they sample regions 
over 100 kpc in linear extent. There are, however, a few clusters containing more 
than one strong, polarized radio source. The cluster Abell 119 (Feretti et al. 
1999) contains three radio galaxies. Using an analysis based on Eq. 3 above, 
Feretti et al. (1999) find that a magnetic field strength of 6-12 u,G extending 
over 3 Mpc could explain the RM distributions for all 3 sources, although they 
note that such a field would exceed the thermal pressure in the outer parts of the 
cluster. 

In a reanalysis of the Abell 119 measurements, Dolag et al. (2001) find that 
the field scales as n^ 9 . This power-law behavior is marginally steeper than that 

2/3 

expected assuming flux conservation, for which the tangled field scales as n e , 
and significantly steeper than that expected assuming a constant ratio between 

1/2 

magnetic and thermal energy density, for which the tangled field scales as n e 
for an isothermal atmosphere. In the 3C 129 cluster there are two extended radio 
galaxies whose RM observations can be fit by a field strength of 6 /iG. Finally in 
A514, Govoni et al. (2001b) has measured the RM distributions of 5 embedded 
(and background) radio sources and found cluster magnetic field strengths of 4-9 
u.G spread over the central 1.4 Mpc of the cluster. If the magnetic field scales 
with the density raised to a positive power, then the product of B and n e in 
Eq. 1 implies that the observed rotation measures are heavily weighted by the 
innermost cells in the cluster (Dreher, Carilli & Perley 1987). 

It has been suggested that high RMs may result from an interaction between 
the radio galaxy and the ICM, such that the RMs are generated locally, and are 
not indicative of cluster magnetic fields. Bicknell, Cameron & Gingold (1990) 
proposed a model in which the RM screen is due to a boundary layer surround- 
ing the radio source in which the large magnetic fields within the radio source 
are mixed with the large thermal densities outside the radio source by Kelvin- 
Helmholtz waves along the contact discontinuity. This model predicts a Faraday 
depolarized region of a few kpc extent surrounding the radio source, where the 
synchrotron emitting material has mixed with the thermal gas. Such a depolar- 
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ized shell has not been observed to date. 

In general, extreme RMs have been observed in sources of very different mor- 
phologies, from edge-brightened (Fanaroff-Riley Class II (Fanaroff & Riley 1974)), 
to edge darkened (FR I) sources. The models for the hydrodynamic evolution of 
these different classes of sources are thought to be very different, with the FRII 
sources expanding supersonically, while the FRI sources expand sub-sonically 
(Begelman, Blandford, & Rees 1984). This argues that the high RMs are not 
solely a phenomenon arising due to a local interaction between the radio source 
and its environment, but are more likely to be a property of the large scale (i.e., 
cluster) environment. Perhaps the most telling argument against the interac- 
tion model is that RM studies of background radio sources seen through cluster 
atmospheres also indicate /iG cluster magnetic fields (see §3.2). 

While we feel that large RMs for cluster center radio sources most likely arise 
in the large scale cluster atmosphere, we should point out that there are some 
cases in which the radio source does appear to compress the gas and fields in the 
ICM to produce local RM enhancements. For example, there is evidence for an 
RM enhancement due to the bow shock preceding the radio hot spots in Cygnus 
A and 3C 194 (Carilli, Perley, & Dreher 1988; Taylor, Inoue, & Tabara 1992). 
However, even in these cases the implied external (i.e., unperturbed) ICM fields 
are a few fiG. 

3.2 Background sources 

The first successful demonstration of Faraday rotation of the polarized emission 
from background radio sources seen through a cluster atmosphere was presented 
by Vallee, MacLeod, & Broten (1986) for A2319. Vallee, MacLeod, & Broten 
(1987) combined the RM excess in A2319 with density estimates from X-ray 
observations by Jones et al. (1979) to estimate a field strength of 2 ^G if the field 
is organized in cells of size 20 kpc. Kim et al. (1990) found that the RMs toward 
sources within 20' of the Coma cluster center had an enhanced RM dispersion by 
38 ± 6 rad m -2 . From this excess they derived a magnetic field strength of 2.0 
± 1.1 /iG assuming a cell size in the range 7-26 kpc. Feretti et al. (1995) found 
evidence from the RM distribution of the embedded cluster source NGC 4869 
for smaller cell sizes (~1 kpc), and subsequently estimated the field strength in 
Coma to be 6.2 /iG. 

The most significant work in this area is the recent VLA survey by Clarke et 
al. (2001), in which they observed radio sources in and behind a representative 
sample of 16 Abell clusters at z < 0.1. They find enhanced rotation measures 
on the large majority of the lines of sight within 0.5 Mpc of the cluster centers 
(Fig. 6), from which they derive an areal filling of the magnetic fields of 95%. 
Their modeling results in magnetic fields of ~5 fiG, assuming a cell size of 10 
kpc. These clusters were chosen not to have cooling flows, but are otherwise 
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Figure 6: The integrated RM plotted as a function of source impact parameter in kiloparsecs 
for the sample of 16 Abell clusters described in Clarke et al (2001). The open symbols represent 
sources viewed through the cluster, while the closed symbols represent the control sample of 
field sources. 



unremarkable in their properties. By observing sources behind the clusters, these 
observations demonstrate that an embedded powerful radio galaxy is not required 
to produce significant RMs. Another advantage of this technique is that it permits 
estimation of the spatial extent of the magnetic fields within the cluster (~ 0.5 
Mpc). The areal filling factor of 95% (assuming constant magnetic fields in and 
among all clusters) suggests a relatively large volume filling factor for the fields, 
with a formal (extreme) lower limit being about 8% for 10 kpc cell sizes. 

3.3 High redshift sources 

Radio galaxies and radio loud quasars have been detected to z = 5.2 (van Breugel 
2000). The extended polarized emission from these sources provides an ideal 
probe of their environments through Faraday rotation observations. Extensive 
radio imaging surveys of z > 2 radio galaxies and quasars have shown large 
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rotation measures, and Faraday depolarization, in at least 30% of the sources, 
indicating that the sources are situated behind dense Faraday screens of mag- 
netized, ionized plasma (Chambers, Miley, & van Breugel 1990: Garrington et 
al. 1988; Carilli, Owen, & Harris 1994; Carilli et al. 1997; Pentericci et al. 2000; 
Lonsdale, Barthel, & Miley 1993; Athreya et al. 1998), with a possible increase 
in this fraction with redshift (Pentericci et al. 2000). Drawing the analogy to 
lower z radio galaxies, these authors proposed that the high z sources may be 
embedded in magnetized (proto-) cluster atmospheres, with \iQ> field strengths. 

A difficulty with the study of high redshift sources is that the sources are 
typically small (< few arcseconds), requiring higher frequency observations (5 
to 8 GHz) in order to properly resolve the source structure. This leads to two 
problems. First is that the rest frame frequencies are then > 20 GHz, such that 
only extreme values of Faraday rotation can be measured (RM > 1000 rad m -2 ). 
And second is that only the flatter spectrum, higher surface brightness radio 
emitting structures in the sources are detected, thereby allowing for only a few 
lines-of-site through the ICM as RM probes. Imaging at frequencies of 1.4 GHz or 
lower with sub-arcsecond resolution is required to address this interesting issue. 

4 Inverse Compton X-ray emission 

Cosmic magnetic fields can be derived by comparing inverse Compton X-ray emis- 
sion and radio synchrotron radiation (Harris & Grindlay 1979; Rephaeli, Gruber, 
& Rothschild 1987). Inverse Compton (IC) emission is the relativisitic extrapola- 
tion of the Sunyaev-ZePdovich effect (Rephaeli 1995), involving up-scattering of 
the ambient photon field by the relativisitic particle population. The IC process 
involves two Lorentz transforms (to and from the rest frame of the electron), plus 
Thompson scattering in the rest frame of the electron, leading to v\q ~ |7 2 ^ g , 
where v\q is the emergent frequency of the scattered radiation, 7 is the elec- 
tron Lorentz factor, and i7, g is the incident photon frequency (Bagchi, Pislar, 
& Lima Neto 1998). From a quantum mechanical perspective, synchrotron ra- 
diation is directly analogous to IC emission, with synchrotron radiation being 
the up-scattering of the virtual photons that constitute the static magnetic field. 
Given a relativistic electron population, the IC emissivity is directly proportional 
to the energy density in the photon field, t/b g , while the synchrotron emissivity is 
proportional to the energy density in the magnetic field, Ub, leading to a simple 
proportionality between synchrotron and IC luminosity: ^7 oc j^-. Given that 
they originate from the same (assumed power-law) relativistic electron popula- 
tion, IC X-rays and synchrotron radio emission share the same spectral index, 
q. The spectral index relates to the index for the power-law electron energy 
distribution, r, as V = 2a — 1, and to the photon index as q — 1. 

In most astrophysical circumstances Ubg is dominated by the cosmic microwave 
background (CMB), except in the immediate vicinity of active star forming re- 
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gions and AGN (Brunetti et al. 2001a; Carilli et al. 2001). The Planck function 
at T = 2.73 K peaks near a frequency of ~ 1.6 x 10 11 Hz, hence IC X- 
rays observed at 20 keV {i>\c = 4.8 x 10 18 Hz), are emitted predominantly by 
electrons at 7 ~ 5000, independent of redshift 1 . The corresponding radio syn- 
chrotron emission from 7 = 5000 electrons peaks at a (rest frame) frequency of 
^ syn - 4.2(j^)7 2 Hz = 100 MHz (Bagchi, Pislar, & Lima Neto 1998). 

Many autnors have considered the problem of deriving magnetic fields by com- 
paring synchrotron radio and inverse Compton X-ray emission (Blumenthal & 
Gould 1970; Harris & Grindlay 1979; Rephaeli, Gruber, & Rothschild 1987; 
Feigelson et al. 1995; Bagchi, Pislar, & Lima Neto 1998). Assuming a = — 1, 
the magnetic field is given by: 

B = 1.7(l + Z ) 2 (f^) - 5 M G (4) 

where S T and S x are the radio and X-ray flux densities at observed frequencies 
v T and respectively. Note that, unlike Faraday rotation measurements, the 
geometry of the field does not play a critical role in this calculation, except in 
the context of the electron pitch angle distribution (see §4.2). 

The principle difficulty in studying IC emission from clusters of galaxies is con- 
fusion due to the thermal emission from the cluster atmosphere. One means of 
separating the two emission mechanisms is through spectroscopic X-ray observa- 
tions at high energy. The IC emission will have a harder, power-law spectrum 
relative to thermal brehmstrahlung emission. Recent high energy X-ray satel- 
lites such as Beppo/Sax and RXTE have revolutionized this field by allowing for 
sensitive observations to be made at energies well above 10 keV (Rephaeli 2001). 
Prior to these instruments, most studies of IC emission from clusters with radio 
halos only provided lower limits to the magnetic fields of about 0.1 u.G (Rephaeli, 
Gruber, & Rothschild 1987). 

Recent observations of four clusters with radio halos with Beppo/Sax and 
RXTE have revealed hard X-ray tails which dominate the integrated emission 
above 20 keV (Rephaeli, Gruber, & Blanco 1999; Fusco-Femiano et al. 2001, 
2000). In Fig. 7 we reproduce the RXTE observations of the Coma cluster. The 
hard X-ray emission in these sources has a spectral index a = — 1.3 ±0.3, roughly 
consistent with the radio spectral index. Comparing the IC X-ray and radio syn- 
chrotron emission in these sources leads to a volume averaged cluster magnetic 
field of 0.2 to 0.4 /j,G, with a relativistic electron energy density ~ 10~ 13 erg 
cm -3 . 

Spatially resolving X-ray observations can also be used to separate non-thermal 
and thermal X-ray emission in clusters. This technique has been used recently 
in the study of the steep spectrum radio relic source in Abell 85 (Bagchi, Pislar, 
& Lima Neto 1998). An X-ray excess relative to that expected for the cluster 

*7 is independent of redshift since !/ bg increases as 1 + z. 
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Figure 7: RXTE spectrum of the Coma cluster. Data and folded Raymond-Smith (kT ~ 7.51 
keV), and power-law (photon index = 2.34) models are shown in the upper frame; the latter 
component is also shown separately in the lower line. Residuals of the fit are shown in the lower 
frame (Rephaeli 2001). 

atmosphere is seen with the ROSAT PSPC detector at the position of the dif- 
fuse radio source in Abell 85 (see Fig. 8). Bagchi, Pislar, & Lima Neto (1998) 
subtract a model of the thermal cluster X-ray emission in order to derive the IC 
contribution, from which they derive a magnetic field of 1.0 ± 0.1 fiG. 

Emission above that expected from the hot cluster atmosphere has also been 
detected in the extreme ultraviolet (EUV = 0.1 to 0.4 keV) in a few clusters 
(Berghofer, Bowyer, & Korpela 2000; Bowyer, Korpela, & Berghofer 2001; Bona- 
mente et al. 2001). It has been suggested that this emission may also be IC 
in origin, corresponding to relativistic electrons with 7 ~ 400 (Atoyan & Volk 
2000). However, the emission spectrum is steep (a < —2), and the EUV emitting 
regions are less extended than the radio regions. Neither of these properties are 
consistent with a simple extrapolation of the radio halo properties to low fre- 
quency (Bowyer, Korpela, & Berghofer 2001). Also, the pressure in this low 7 
relativistic component would exceed that in the thermal gas by at least a factor 
of three (Bonamente et al. 2001). 
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R.A. 

Figure 8: The cluster Abell 85 central region at different wavelengths. The photographic R- 
filter image (UK Schmidt Telescope and the Digitized Sky Survey) is shown in grey scale. Pull 
contour lines show the multiscale wavelet reconstruction of the ROSAT PSPC X-ray data shown 
in Fig. 1. The OSRT 326.5 MHz radio surface- brightness image is depicted using dot-dashed 
contour lines. All contours are spaced logarithmically (Bagchi, Pislar, & Lima Neto 1998). 
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4-1 Electron lifetimes 

An important point concerning IC and synchrotron emission from clusters is 
that of particle lifetimes. The lifetime of a relativistic electron is limited by IC 
losses off the microwave background to: tic = 7.7 x 10 9 (^)(1 + z)~ 4 years 
(Sarazin 2001a). 2 Relativistic electrons emitting in the hard X-ray band via 
IC scattering of the CMB have lifetimes of about 10 9 years, while the lifetimes 
for 1.4 GHz synchrotron emitting electrons are a factor of four or so shorter in 
fiG fields. Diffusion timescales (set by streaming along the tangled magnetic 
field) for cluster relativistic electrons are thought to be longer than the Hubble 
time (Sarazin 2001a; Colafrancesco & Blasi 1998; Casse, Lemoine & Pelletier 
2001), making cluster atmospheres efficient traps of relativistic electrons, much 
like galaxy disks. The fact that the diffusion timescales are much longer than the 
energy loss timescales for 7 > 10 4 electrons requires in situ acceleration in order 
to explain radio halo sources (Schlikeiser, Sievers, & Thiemann 1987). 

Cluster merger shock fronts are obvious sites for first order Fermi accelera- 
tion, while subsequent turbulence may lead to second order (stochastic) Fermi 
acceleration (Brunetti et al. 2001b; Eilek 1999; Ensslin et al. 1998; Markevitch 
& Vikhlinin 2001). Active particle acceleration during cluster mergers provides 
a natural explanation for the observed correlation between cluster radio halos 
and substructure in cluster atmospheres (Govoni et al. 2001a), and between clus- 
ter radio luminosity and cluster atmosphere temperature, assuming that the gas 
temperature increases during a merger (Liang et al. 2000). Brunetti et al. (2001a) 
develop a two phase model in which initial relativistic particle injection into the 
ICM occurs early in the cluster lifetime by starburst driven winds from cluster 
galaxies, and/or by shocks in early sub-cluster mergers. The second phase in- 
volves re-acceleration of the radiatively aged particle population via more recent 
cluster mergers. Their detailed application of this model to the Coma cluster 
suggests a merger has occurred within the last 10 9 years. 

Another mechanism proposed for in situ relativistic particle injection is sec- 
ondary electron production via the decay of 7r-mesons generated in collisions 
between cosmic ray ions (mostly protons) and the thermal ICM (Dennison 1980; 
Ensslin & Gopal-Krishna 2001). The important point in this case is that the ener- 
getic protons have radiative lifetimes orders of magnitude longer than the lower 
mass electrons. The problem with this hypothesis is that the predicted 7-ray 
fluxes exceed limits set by EGRET by a factor of 2 to 7 (Blasi & Colafrancesco 
1999). 

2 For 7 > 300, IC losses dominate (or synchrotron losses for B > 2.3(1 + z) 2 /<G), while for 
lower 7 electrons Brehmstrahlung losses dominate in cluster environments (Sarazin 2001b). 
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4-2 Reconciling IC and RM- derived fields 

The IC estimated cluster magnetic fields are typically 0.2 to 1 /iG, while those ob- 
tained using RM observations are an order of magnitude higher. Petrosian (2001) 
has considered this discrepancy in detail, and finds that the different magnetic 
field estimates can be reconciled in most cases by making a few reasonable as- 
sumptions concerning the electron energy spectrum and pitch angle distribution. 

First, an anisotropic pitch angle distribution biased toward low angles would 
clearly weaken the radio synchrotron radiation relative to the IC X-ray emission. 
Such a distribution will occur naturally due to the fact that electrons at large pitch 
angles have greater synchrotron losses. A potential problem with this argument 
is that pitch-angle scattering of the relativistic electrons by Alfven waves self- 
induced by particles streaming along field lines is thought to be an efficient process 
in the ISM and ICM (Wentzel 1974), such that re-isotropization of the particle 
distribution will occur on a timescale short compared to radiative timescales. 
Petrosian points out that most derivations of magnetic fields from IC emission 
assume the electrons are gyrating perpendicular to the magnetic field lines. Just 
assuming an isotropic relativistic electron pitch angle distribution raises the IC- 
estimated magnetic field by a factor of two or so. 

And second, the IC hard X-ray emission is from relativistic electrons with 
7 ~ 5000, corresponding to radio continuum emission at 100 MHz for /iG mag- 
netic fields. Most surveys for cluster radio halos have been done at 1.4 GHz 
(Giovannini, Tordi & Feretti 1999; Govoni et al. 2001a), corresponding to elec- 
tron Lorentz factors of 7 ~ 18000. A steepening in the electron energy spectrum 
at Lorentz factors around 10 4 will reduce the 1.4 GHz radio luminosities, but 
retain the IC hard X-ray emission. For example, Petrosian finds that a steepen- 
ing in the power-law index for the particle energy distribution from T = — 3 to 
r = —5 (corresponding to a = —1 to —2) at 7 ~ 10 raises the IC-estimated 
fields to ~ 1/iG. Such a steepening of the electron energy spectrum at 7 ~ 10 4 
will arise naturally if no relativistic particle injection occurs over a timescale 
~ 10 8 years (see §4.1). The problem in this case is the fine tuning required to 
achieve the break in the relevant energy range for a number of clusters. In gen- 
eral, a negatively curved (in log-space) electron energy distribution will inevitably 
lead to IC estimated fields being lower than those estimated from 1.4 GHz radio 
observations, unless a correction is made for the spectral curvature. 

Others have pointed out that magnetic substructure, or filamentation, can lead 
to a significant difference between fields estimated using the different techniques. 
A large relativistic electron population can be 'hidden' from radio observations 
by putting them in weak field regions (Feretti et al. 2001; Rephaeli, Gruber, & 
Rothschild 1987; Rephaeli 2001; Goldshmidt & Rephaeli 1993; Rudnick 2000). 
A simple example of this is if the relativistic particles have a larger radial scale- 
length than the magnetic fields in the cluster. In this case, most of the IC 
emission will come from the weak field regions in the outer parts of the cluster, 
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while most of the Faraday rotation and synchrotron emission occurs in the strong 
fields regions in the inner parts of the cluster. 

Another explanation for the discrepancy between IC and RM derived magnetic 
fields is to assume that the hard X-rays are not IC in origin, in which case the 
IC estimates become lower limits. A number of authors have considered high 
energy X-ray emission due to non-thermal Brehmstrahlung, i.e. Brehmstrahlung 
radiation from a supra-thermal tail of electrons arising via stochastic acceleration 
in a turbulent medium (Fermi II acceleration) (Blasi 2000; Sarazin 2001a; Ensslin 
& Gopal-Krishna 2001; Dogiel 2000). The problem with this hypothesis is the 
energetics: Brehmstrahlung is an inefficient radiation mechanism, with most of 
the collisional energy going into heat. Assuming an energy source is available 
to maintain the supra-thermal tail, Petrosian (2001) shows that the collisional 
energy input by the supra-thermal particles would be adequate to evaporate the 
cluster atmosphere on a timescale of order 10 8 years. 

The current hard X-ray spectroscopic observations are limited to very low 
spatial resolution (~ 1°), while X-ray imaging instruments have high energy cut- 
offs at around 10 keV. Likewise, sensitive, arcminute resolution radio images for a 
large number of clusters are available only at 1.4 GHz, corresponding to electrons 
with Lorentz factors 3 to 4 times higher than those emitting hard X-rays. Both of 
these limitations will be overcome in the coming years with the launch of hard X- 
ray imaging satellites such as Constellation-X, and improvements in radio imaging 
capabilities at 300 MHz and below at the Very Large Array and the Giant Meter 
Wave Radio Telescope. 

5 Cold fronts 

It has long been noted that in order to maintain temperature gradients in X-ray 
clusters the thermal conduction must be suppressed by two orders of magnitude 
relative to the classical Spitzer value (McKee & Begelman 1990; Fabian 1994; 
Spitzer 1962). If not, cooler structures on scales ~ 0.1 Mpc will be evaporated 
by thermal conduction from the hot surrounding medium on timescales ~ 10 7 
years. Examples of such cooler structures in clusters include cooling flow cluster 
cores and X-ray corona surrounding large galaxies (Fabian 1994; Vikhlinin et al. 
2001b), with temperature differences ranging from a factor of 2 to 5 relative to 
the hot ICM. 

Cowie & McKee (1977) show that the conductivity can be suppressed by almost 
an order of magnitude below the Spitzer value due to the development of electro- 
static fields in cases where the Coulomb mean free path (mfp) is comparable to 
the scale of thermal gradients. For large scale structure in cluster atmospheres 
this reduction is not adequate since the mfp ~ 1.2 x 10 22 ( 5x ^7 K ) 2 ( ooicm- ;i ) _1 
cm, or just a few kpc for a typical cluster. 

The presence of magnetic fields will reduce the conductivity in a thermal plasma 
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(Field 1965; Parker 1979; Binney & Cowie 1981; Chandran, Cowley, & Albright 
1999). The simple point is that the gyro radius for thermal electrons in the ICM 
is ~ 2 x 10 8 ( jj^-r )( 5x ^ 7 K ) cm, many orders of magnitude below the collisional 
mfp. Tribble (1993) shows that the presence of a cluster magnetic field will 
lead naturally to the development of a multiphase ICM, with thermally isolated 
regions on scales set by the magnetic structures (although cf. Rosner & Tucker 
(1989)). 

The idea of magnetic suppression of thermal conductivity in cluster gas has 
been verified with the recent discovery of 'cold fronts' in clusters of galaxies 
(Markevitch et al. 2000; Vikhlinin et al. 2001b, a). These fronts manifest them- 
selves as sharp discontinuities in X-ray surface brightness (Fig. 9). They are not 
shocks, since the increase in density is accompanied by a decrease in tempera- 
ture such that there is no dramatic change in the pressure and entropy across 
the front (Markevitch et al. 2000; Ettori & Fabian 2000). These structures are 
interpreted as resulting from cluster mergers, where a cooler subcluster core falls 
into a hot ICM at sub- or trans-sonic velocities (~ W fi km s _1 ). A discontinuity 
is formed where the internal pressure of the core equals the combined ram and 
thermal pressure of the external medium. Gas beyond this radius is stripped 
from the merging subcluster, and the core is not penetrated by shocks due to its 
high internal pressure. 

The best example of a cluster cold front is that seen in Abell 3667 (see Fig. 8) 
(Vikhlinin et al. 2001b, a). In this case, the temperature discontinuity occurs over 
a scale of ~ 5 kpc, comparable to the collisional mfp, thereby requiring thermal 
isolation. Magnetic fields play a fundamental role in allowing for such structures 
in two ways: (i) by suppressing thermal conduction, and (ii) by suppressing 
Kelvin-Helmholtz mixing along the contact discontinuity. Vikhlinin et al. (2001a) 
present a model in which the field is tangentially sheared by fluid motions along 
the contact discontinuity. They invoke magnetic tension to suppress mixing, and 
show that the required magnetic pressure is between 10% and 20% of the thermal 
pressure. The implied fields are between 7 and 16 u.G. They also argue that the 
fields cannot be much stronger than this, since dynamically dominant fields would 
suppress mixing along the entire front, which does not appear to be the case. 

The existence of cold fronts provides strong evidence for cluster magnetic fields. 
However, the field strengths derived correspond to those in the tangentially 
sheared boundary region around the front. Relating these back to the unper- 
turbed cluster field probably requires a factor of a few reduction in field strength, 
implying unperturbed field strengths between 1 and 10 fiG, although the exact 
scale factor remains uncertain (Vikhlinin et al. 2001a). 
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Figure 9: (a) A smoothed 0.5-4 keV Chandra image of Abell 3667. The most prominent 
feature is the sharp surface brightness edge (cold front). The front shape is nearly circular as 
indicated by the arc. (b) Temperature map. The typical statistical error in this image is ±1 
keV. The cold, ~ 4 keV, region near the center of the map coincides with the inside of the cold 
front (Vikhlinin et al. 2001a). 



24 



Cluster Magnetic Fields 



6 Magnetic support of cluster gas and the baryon crisis 

Two interesting issues have arisen in the study of cluster gas and gravitational 
masses. First is the fact that total gravitating masses derived from weak grav- 
itational lensing are a factor of a few higher than those derived from X-ray ob- 
servations of cluster atmospheres assuming hydrostatic equilibrium of isothermal 
atmospheres (Loeb & Mao 1994; Miralda-Escude & Babul 1995). And second is 
the 'baryon crisis', in which the baryonic mass of a cluster, which is dominated by 
the mass of gas in the hot cluster atmosphere, corresponds to roughly 5% of the 
gravitational mass derived assuming hydrostatic equilibrium for an isothermal 
cluster atmosphere. This baryon fraction is a factor of three to five larger than 
the baryon fraction dictated by big bang nucleosynthesis in inflationary world 
models (White et al. 1993). 

A possible solution to both these problems is to invoke non-thermal pressure 
support for the cluster atmosphere, thereby allowing for larger gravitating masses 
relative to those derived assuming hydrostatic equilibrium. A number of authors 
have investigated the possibility of magnetic pressure support for cluster atmo- 
spheres (Loeb & Mao 1994; Miralda-Escude & Babul 1995; Dolag & Schindler 
2000). The required fields are about 50^G, which is an order of magnitude, or 
more, larger than the measured fields in most cases, except perhaps in the inner 
10's of kpc of cooling flow clusters. For most relaxed clusters Dolag & Schindler 
(2000) find that magnetic pressure affects hydrostatic mass estimates by at most 
15%. 

Other mechanisms for non-thermal pressure support of cluster atmospheres in- 
volve motions of the cluster gas other than thermal, such as turbulent or bulk 
motions due to a recent cluster merger (Mazzotta et al. 2001; Wu 2000). For 
relaxed clusters a number of groups have shown that the lensing and X-ray mass 
estimates can be reconciled by using non-isothermal models for the cluster at- 
mospheres, i.e., by allowing for radial temperature gradients (Allen, Ettori, & 
Fabian 2001b; Markevitch et al. 1999). 

7 GZK limit 

Greisen (1966) and Zatsepin & Kuzmin (1966) pointed out that EeV cosmic rays 
lose energy due to photo-pion production through interaction with the CMB. 
These losses limit the propagation distance for 10 20 eV particles to about 60 
Mpc. Yet no clear correlation has been found between the arrival direction of 
high energy cosmic rays and the most likely sites of origin for EeV particles, 
namely, AGN at distances less than 60 Mpc (Elbert & Sommers 1995; Biermann 
1999). One solution to the GZK paradox is to assume the energetic particles 
are isotropized in the IGM by tangled magnetic fields, effectively randomizing 
their observed arrival direction. Such isotropization requires fields in the local 



Carilli & Taylor 



25 



super-cluster > 0.3/xG (Farrar 2000; Isola, Lemoine, & Sigl 2001). 
8 Synthesis 

In Table 1 we summarize the cluster magnetic field measurements. Given the 
limitations of the current instrumentation, the limited number of sources studied 
thus far, and the myriad physical assumptions involved with each method, we are 
encouraged by the order-of-magnitude agreement between cluster field strengths 
derived from these different methods. Overall, the data are consistent with cluster 
atmospheres containing ~ u.G fields, with perhaps an order of magnitude scatter 
in field strength between clusters, or within a given cluster. 

The rotation measure observations of background radio sources, and in partic- 
ular the observations of a complete X-ray selected sample of clusters by Clarke, 
Kronberg, & Bohringer (2001) dictate that jj,G magnetic fields with high areal 
filling factors are a standard feature in clusters, and that the fields extend to 
large radii (0.5 Mpc or more). The rotation measure observations of extended 
radio galaxies embedded in clusters impose order on the fields, with coherence 
scales of order 10 kpc, although larger scale coherence in overall RM sign can 
be seen in some sources. Observations of inverse Compton emission from a few 
clusters with radio halos provide evidence against much stronger, pervasive and 
highly tangled fields. 

In most clusters the fields are not dynamically important, with magnetic pres- 
sures one to two orders of magnitude below thermal gas pressures. But the fields 
play a fundamental role in the energy budget of the ICM through their effect on 
heat conduction, as is dramatically evident in high resolution X-ray observations 
of cluster cold fronts. 

If most clusters contain u.G magnetic fields, then why don't most clusters have 
radio halos? The answer may be the short lifetimes of the relativistic electrons 
responsible for the synchrotron radio emission (see §4.1). Without re-acceleration 
or injection of relativistic electrons, a synchrotron halo emitting at 1.4 GHz will 
fade in about 10 years due to synchrotron and inverse Compton losses. This may 
explain the correlation between radio halos and cluster mergers, and the anti- 
correlation between radio halos and clusters with relaxed X-ray morphologies. In 
this case, the fraction of clusters with radio halos should increase with decreasing 
survey frequency. 

The existence of /iG-level fields in cluster atmospheres appears well established. 
The challenge for observers now becomes one of determining the detailed prop- 
erties of the fields, and how they relate to other cluster properties. Are the fields 
filamentary, and to what extent do the thermal and non-thermal plasma mix in 
cluster atmospheres? What is the radial dependence of the field strength? How 
do the fields depend on cluster atmospheric parameters, such as gas temperature, 
metalicity, mass, substructure, or density profile? How do the fields evolve with 
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Table 1: Cluster Magnetic Fields 



Method 


Strength 


Model Parameters 








Synchrotron Halos 


0.4 - 1 


minimum energy, k = 7} — 1, 






i^iow = 10MHz, v high = 10GHz 


Faraday rotation (embedded) 


3-40 


cell size = 10 kpc 


Faraday rotation (background) 


1 - 10 


cell size — 10 kpc 


Inverse Compton 


0.2 - 1 


a = -1, 7radio ~ 18000, 7xr»y ~ 5000 


Cold Fronts 


1 - 10 


amplification factor ~ 3 


GZK 


> 0.3 


AGN = site of origin for EeV CRs 



cosmic time? And do the fields extend to even larger radii, perhaps filling the 
IGM? The challenge to the theorists is simpler: how were these field generated? 
This topic is considered briefly in the next section. 

9 Field origin 

When attempting to understand the behavior of cosmic magnetic fields, a critical 
characteristic to keep in mind is their longevity. The Spitzer conductivity (Spitzer 
1962) of the ICM is: a ~ 3x 10 sec -1 (for comparison, the conductivity of liquid 
mercury at room temperature is 10 sec -1 ). The timescale for magnetic diffusion 
in the ICM is then: r dif f = 4na(^) 2 ~ W 36 (y$^) 2 years, where L is the spatial 
scale for magnetic fluctuations. The magnetic Reynold's number is: R m = ~ 

1Q29 (iOkl^)( ioookm s-i )' wnere T"conv = the convective timescale = ^, and V is the 
bulk fluid velocity. The essentially infinite diffusion timescale for the fields implies 
that once a field is generated within the ICM, it will remain extant unless some 
anomalous resistive process occurs e.g., reconnection via plasma wave generation 
in shocks. 

Perhaps the simplest origin for cluster magnetic fields is compression of an 
intergalactic field. Clusters have present day overdensities 5 ~ 10 3 . In order to 
get -Bicm > 10 _7 G by adiabatic compression (B oc 5~s) then requires IGM fields 
Sigm > 10"VG. 

Of course, this solution merely pushes the field origin problem from the ICM 
into the IGM. An upper limit to IGM fields of 10 _9 G is set by Faraday rotation 
measurements of high z radio loud QSOs, assuming a cell size of order 1 Mpc 
(Kronberg 1996; Blasi, Buries, & Olinto 1999). A limit to IGM magnetic fields at 
the time of recombination can also be set by considering their affect on the CMB. 
Dynamically significant magnetic fields will exert an anisotropic pressure on the 
gas, which must be balanced by gravity. Detailed studies of this phenomenon 
in the context of recent measurements of the CMB anisotropies shows that the 
comoving IGM fields 3 must be less than a fewxlO -9 G (Barrow, Ferreira, & Silk 

3 Comoving fields correspond to equivalent present epoch field strengths, i.e. corrected for 
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1997; Clarkson & Coley 2001; Adams et al. 1996). A co-moving field of 10" 9 G 
at recombination would lead to Faraday rotation of the polarized CMB emission 
by 1° at an observing frequency of 30 GHz, a measurement that is within reach 
of future instrumentation (Kosowsky & Loeb 1996; Grasso and Rubinstein 2001). 
Considerations of primordial nucleosynthesis and the affect of magnetic fields on 
weak interactions and electron densities imply upper limits to comoving IGM 
fields of 10 -7 G (Grasso and Rubinstein 1995). 

The origin of IGM magnetic fields has been considered by many authors. One 
class of models involves large scale field generation prior to recombination. An 
excellent review of pre-recombination magnetic field generation is presented by 
(Grasso and Rubinstein 2001). Early models for pre-recombination field genera- 
tion involved the hydrodynamical ('Biermann') battery effect (Biermann 1950). 
In general, the hydrodynamic battery involves charge separation arising from the 
fact that electrons and protons have the same charge, but very different masses. 
For instance, protons will have larger Coulomb stopping distances than electrons, 
and be less affected by photon drag. Harrison (1970) suggested that photon drag 
on protons relative to electrons in vortical turbulence during the radiation era 
could lead to charge separation, and hence magnetic field generation by electric 
currents. Subsequent authors have argued strongly against vortical density per- 
turbations just prior to recombination, since vortical (and higher order) density 
perturbations decay rapidly with the expansion of the universe (Rees 1987). This 
idea has been revisited recently in the context of vortical turbulence generated by 
moving cosmic strings (Vachaspati & Vilenkin 1991; Avelino & Shellard 1995). 
Other mechanisms for field generation prior to recombination include battery 
affects during the quark-hadron (QCD) phase transition (Quashnock, Loeb, & 
Spergel 1989), dynamo mechanisms during the electro- weak (QED) phase tran- 
sition (Baym, Bodeker, & McLerran 1996), and mechanisms relating to the basic 
physics of inflation (Turner & Widrow 1988). 

A problem with all these mechanisms is the survivability of the fields on relevant 
scales during the radiation era. Battaner & Lesch (2000) argue that magnetic 
and photon diffusion will destroy fields on comoving scales < few Mpc during 
this epoch, thereby requiring generation of the fields in the post-recombination 
universe by normal plasma processes during proto-galactic evolution (see also 
Lesch & Birk (1998)). 

Models for post-recombination IGM magnetic field generation typically involve 
ejection of the fields from normal or active galaxies (Kronberg 1996; Rees 1989). 
A simple but cogent argument in this case is that the metalicity of the ICM is 
typically about 30% solar, implying that cluster atmospheres have been polluted 
by outflows from galaxies (Aguirre et al. 2001). A natural extension of this 
idea would be to include magnetic fields in the outflows (Goldshmidt & Rephaeli 
1993). It has also been suggested that IGM fields could be generated through 

cosmic expansion assuming flux freezing. 
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turbulent dynamo processes and/or shocks occurring during structure formation 
(Zweibel 1988; Kulsrud et al. 1997; Waxman & Loeb 2000), or by battery effects 
during the epoch of reionization (Gnedin et al. 2000). 

Seed magnetic fields will arise in the earliest stars via the normal gas kinemati- 
cal Biermann battery mechanism. These fields are amplified by the a — Q dynamo 
operating in stellar convective atmospheres (Parker 1979), and then are ejected 
into the ISM by stellar outflows and supernova explosions. The ISM fields can 
then be injected into the IGM by winds from active star forming galaxies (Heck- 
man 2001). Kronberg, Lesch, & Hopp (1999) consider this problem in detail and 
show that a population of dwarf starburst galaxies at z > 6 could magnetize 
almost 50% of the universe, but that at lower redshifts the IGM volume is too 
large for galaxy outflows to affect a significant fraction of the volume. 

De Young (1992) and Rephaeli (1988) show that galaxy outflows, and/or gas 
stripping by the ICM, in present day clusters are insufficient to be solely respon- 
sible for cluster fields ~ 1/xG without invoking subsequent dynamo amplification 
of the fields by about an order of magnitude in the cluster atmosphere. A simple 
argument in this case is that the mean density ratio of the ICM versus the ISM, 
S ~ 0.01, such that ICM fields would be weaker than ISM fields by <5§ ~ 0.05, 
corresponding to maximum ICM fields of 0.2 to 0.5 fiG. 

Fields can be ejected from Active Galactic Nuclei (AGN) by relativistic out- 
flows (radio jets) and Broad Absorption Line outflows (BALs) (Rees & Setti 1968; 
Daly & Loeb 1990). The ultimate origin of the fields in this case may be a seed 
field generated by a gas kinematic battery operating in the dense accretion disk 
around the massive black hole, plus subsequent amplification by an a — SI dynamo 
in the rotating disk (Colgate & Li 2000). Detailed consideration of this problem 
(Furlanetto & Loeb 2001; Kronberg et al. 2001) using the statistics for high z 
QSO populations shows that by z ~ 3, between 5% and 20% of the IGM may be 
permeated by fields with energy densities corresponding to > 10% the thermal 
energy density of the photo-ionized IGM at 10 4 K, corresponding to comoving 
field strengths of order 10 -9 /iC 

Kronberg et al. (2001) point out that powerful double radio sources such as 
Cygnus A (radio luminosities ~ 10 45 erg s -1 ) typically have total magnetic ener- 
gies of about 10% that of the ICM as a whole. Hence, about ten powerful double 
radio sources over a cluster lifetime would be adequate to magnetize the cluster 
at the /iG level. 

Galaxy turbulent wakes have been proposed as a means of amplifying cluster 
magnetic fields (Jaffe 1980; Tribble 1993; Ruzmaikin, Sokolov, & Shukurov 1989). 
The problem in this case is that the energy appears to be insufficient, with ex- 
pected field strengths of at most ~ 0.1/iG. Also, the size scale of the dominant 
magnetic structures is predicted to be significantly smaller than the 5 to 10 kpc 
scale sizes observed (Goldshmidt & Rephaeli 1993; De Young 1992). 

Cluster mergers are the most energetic events in the universe since the big bang, 
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releasing of order 10 ergs in gravitational binding energy (Sarazin 2001a). For 
comparison, the total thermal energy in the cluster atmosphere is 
~ 10 63 ( io-um, X sxio^k ) er S s ; an d the total energy contained in the cluster mag- 
netic fields is ~ 10 60 (j^) 2 ergs. Hence, only a fraction of a percent of the 
cluster merger energy need be converted into magnetic fields. One possibility 
for merger-generated magnetic fields is a rotational dynamo associated with he- 
lical turbulence driven by off-center cluster mergers. This mechanism requires 
net cluster rotation - a phenomenon that has yet to be seen in cluster galaxy 
velocity fields (cf. (Dupke & Bregman 2001)). The lack of observed rotation 
for clusters suggests low impact parameters for mergers (< 100 kpc) on average 
(Sarazin 2001a), as might arise if most mergers occur along filamentary large 
scale structure (Evrard & Gioia 2001). The energetics of even slightly off-center 
cluster mergers is adequate to generate magnetic fields at the level observed, but 
the slow cluster rotation velocities (< 100 km s _1 ) imply only one or two rota- 
tions in a Hubble time (Colgate & Li 2000), which is insufficient for mean field 
generation via the inverse cascade of the a — Q dynamo (Parker 1979). 

A general treatment of the problem of magnetic field evolution during cluster 
formation comes from numerical studies of heirarchical merging of large scale 
structure including an initial intergalactic field ~ 10 -9 G ((Dolag & Schindler 
2000; Roettiger, Stone, & Burns 1999). These studies show that a combination 
of adiabatic compression and non-linear amplification in shocks during cluster 
mergers may lead to ICM mean fields of order lfiG. 

A related phenomenon is field amplification by (possible) cooling flows. Soker 
& Sarazin (1990) have considered this mechanism in detail, and show that the 
amplification could be a factor of ten or larger in the inner 10's of kpc. They 
predict a strong increase in RMs with radius (oc r 2 ), with centrally peaked radio 
halos. Such an increase may explain the extreme RM values seen in powerful radio 
sources at the centers of cooling flow clusters (see §3.1), although the existence 
of gas inflow in these systems remains a topic of debate (Binney 2001). 

Overall, there are a number of plausible methods for generating cluster mag- 
netic fields, ranging from injection of fields into the IGM (or early ICM) by active 
star forming galaxies and/or radio jets at high redshift, to field amplification by 
cluster mergers. It is likely that a combination of these phenomena give rise to 
the /iG fields observed in nearby cluster atmospheres. Tests of these mechanisms 
will require observations of (proto-) cluster atmospheres at high redshift, and a 
better understanding of the general IGM field. 
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